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ABSTRACT 

We study the dependence of quasar clustering on quasar luminosity and black hole mass by mea¬ 
suring the angular over density of photometrically selected galaxies imaged by WISE about z ~ 0.8 
quasars from SDSS. By measuring the quasar-galaxy cross-correlation function and using photomet¬ 
rically selected galaxies, we achieve a higher density of tracer objects and a more sensitive detection 
of clustering than measurements of the quasar autocorrelation function. We test models of quasar 
formation and evolution by measuring the luminosity dependence of clustering amplitude. We find 
a significant overdensity of WISE galaxies about z ^ 0.8 quasars at 0.2-6.4 h~^ Mpc in projected 
comoving separation. We find no appreciable increase in clustering amplitude with quasar luminosity 
across a decade in luminosity, and a power-law fit between luminosity and clustering amplitude gives 
an exponent of —0.01 ± 0.06 (1 a errorbar). We also fail to find a significant relationship between 
clustering amplitude and black hole mass, although our dynamic range in true mass is suppressed due 
to the large uncertainties in virial black hole mass estimates. Our results indicate that a small range 
in host dark matter halo mass maps to a large range in quasar luminosity. 

Subject headings: quasars: general, galaxies: active, large-scale structure of universe 


1. BACKGROUND 


Direct observations have recently confirmed that 
nearly every nearby elliptical galaxy or spiral bulge hosts 
a central s upermassive black hole (Kormendy &: Rich- 
stone|[T9^ I. Further studies have shown that black hole 


properties are tightly related to the macroscopic proper¬ 
ties of the host galaxy. Black hole mass is closely tied to 


both the stellar velocity dispersion (Gebhardt et al. 2000 


k HuntpboS 

Magorrian et al.||1998l 

. Both the velocity 

dispersion anc 

spheroid mass are rel 

ited to the 

mass of 

the host dark matter halo (Ferrarese 

20021. These rela- 

tionships provide evidence for a link 

between black hole 

and host galaxy evolution (Silk k Rees 

1998 

Hopkins 


Local s upermassive black holes are relics of past quasar 
activity (Soltan 19821. Quasars are an important stage 
in the coevolution of a black hole and its host galaxy. 
Quasar activity is associated with massive, disruptive 
flows of gas into the center of the galaxy, and simula¬ 
tions predict that quasar activity suppresses star forma¬ 
tion, leading to a link between black h ole mass and host 
spheroid mass ( Di Matteo et al.]|2005 ). 

However, the details of quasar formation and evolution 
are not fully known. Many different models have been 
proposed to explain quasar activity, ranging in sophis¬ 
tication from numerical sim ulations of quasar activity 
tracking gas hydrodynamics dHopkins et al.|2008 IDegraf 


et al.|2011 Thacker et al.|200y Chatterjee et al.|2U12p to 

semi-analytical methods tracking only dark matter halos 

jKau^m^^&Haehnelt 2000 2002 Bonoli et al. 2009 


20101. A brief, comprehensive, and useful overview of 
various quasar mod els can be found in Appendix B of 
White et al. (2012). In this paper, we consider the dif- 
terences between low and high luminosity quasars, and 
between low and high mass black holes. Specifically, we 
study the relationships between luminosity and host halo 
mass and between black hole mass and host halo mass 
by measuring quasar clustering. 

Quasar clustering is a valuable tool for exploring 
quasar properties and evolution. Quasar clustering is 
often measured using the linear quasar bias bg, defined 
by 

— ^Q^matter (1) 

where is the quasar autocorrelation function and 
^matter is the Correlation function of the linear-regime 
matter field. The bias is then converted to the host halo 
mass usi ng analytic formulas tested against iV-body sim¬ 


ulations (Cole & Kaiser 1989 Mo & White 


et al.|20()T| 7 which predict that at fixed redsnitt the halo 


1996 


Sheth 


mass IS a monotonically increasing function of bias. Mea¬ 
surements of the luminosity dependence of quasar clus¬ 
tering allow estimation of the slope and scatter of the re¬ 
lationship between luminosity and host halo mass. The 
halo mass-luminosity relationship arises from the combi¬ 
nation of a black hole mass-luminosity relationship and 
a halo mass-black hole mass relationship. By also mea¬ 
suring the black hole mass dependence of clustering, we 
can separate the contributions of the black hole mass- 
luminosity relationship from the halo mass-black hole 
mass relationship. 

The quasar lifetime can be determined from measure- 


Eanidakis et al.|2012[ ) to scaling r elations with no p refer- ments of the quasar bias ( ^le&Ka^CT 


ences tor the uriderlying physics (Conroy fc White 2013 
|Lidz et al.|2006[[Croton[2009||bhen|2UU9HBooth fc Schaye 
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& Hui 2001 Martini fc Weinberg| 2001 ~ 


1989 


Haiman 


quasars 

are long-lived, then most of the quasars in the universe 
are actively accreting and quasars are rare, whereas if 
quasars are short-lived, there are many unobserved dor- 
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mant quasars for each observed active quasar and quasars 
are c ommon. Since high ma ss, highly biased halos are 
rare (Press & Schechter|l9741 a highly biased quasar pop- 
ulation implies a long litetiine. The conversion between 
bias and lifetime assumes that luminosity is a monotonic 
function of halo mass with no scatter. If there is scatter 
in the luminosity-halo mass relatio nship, the lifetime is 


et al.||2008[ IShen et al. |2009i IMountrichas et al.||200^ 


great er than for the no-scatter case (Martini & Weinberg 
20011. Thus, measurements of the luminosity depen¬ 
dence of quasar clustering allow for joint determination of 
the quasar lifetime and the scatter in th e luminosity-halo 
mass relationship ( Sh ankar et al. ||2010t. The ana l ysis of 
Martini & Weinberg (|200l| andjjiainian & Hui (2001| 
assumed that the quasar lightcurve is a step function, ei¬ 
ther “on” or “off” at any given time. More sophisticated 
models of the quasar lightcurve have been pro posed (e.g. 
Kauffmann & Haehnelt 2002||Lidz et al.|2006|), leading to 
predictions ot luminosity or redshitt-dependent lifetimes. 
In principle, measurements of the luminosity dependence 
of clustering can help distinguish between different mod¬ 
els of the quasar lightcurve. 

Large surveys su ch as the Sloan Digital Sky Survey 


(York et al. 2000|and the 2dF QSO Redshift Survey 


(Croom et al. '|2UU4 1 provide the large samples neces- 


sary to study the luminosity dependence of quasar clus¬ 
tering. Early studies measured the luminosity depen¬ 
dence of the real-space and redshift-space autocorrela- 
tion function using spectroscopically selected quasars: 
Croom et al. (20051 found no luminosity dependence of 
the redshitt-space clustering strength of 2dF q uasars at 
0.5 < z < 2.5, and Porciani & Norberg (2006| found no 
luminosity dependence on real-space clustering strength 
for 2dF quasars at 0.8 < z < 1.3 and only marginal 
evidence at 1.3 < z < 2.1. More recently, [Shanks] 


et al. (2011) found no luminosity dependence on clus¬ 
tering strength at z = 1.4 using quasar autocorrelation 
measuremen ts from 2QZ, 2SLAQ, and SDSS, and [White] 
et al. (2012) found no luminosity dependence on redshift- 
space clustering strength for SDSS quasars at z = 2.4. 
Due to the low space density of spectroscopically de¬ 
tected quasars, these studies generally suffer from large 
error bars. One way to obtain smaller error bars is to 
measure the auto correlation function of photometrically 
selected quasars. Myers et al. (2007) measured the au¬ 
tocorrelation of 30U,(J(JU photometric quasars from SDSS 
and found no luminosity dependence of quasar bias at 
photometric redshifts z = 0.85, z = 1.44, and z = 1.92. 

Another approach is to measure the cross-correlation 
of galaxies about quasars. Since galaxies have a much 
higher spatial density than quasars, cross-correlations 
can achieve lo wer error bars than autocorrela tions (as 
suggested by Kauffmann & Haehnelt 20021. Early 
quasar-galaxy cross-correl ations used small samp l es and 
thus had large error bars. Adelberger & Steidelf (20051, 
measuring the cross-correlation ot quasars at 2 < z < 3 
with Lyman-break galaxies, found no luminosity depen¬ 
dence across 4 decades in quasar luminosity, but urged 
caution in inter pretation of their result due to large error 
bars. Similarly, Coil et al. (2007) also found no luminos¬ 
ity dependence on quasar-galaxy clustering at z « 1, but 
again with relatively large error bars. More recent studies 
found little evidence for luminosity-dependen t clustering 


Shirasaki et al. 2011 Shen et al. 2013p, although Shen 

et al. ~( ^009| did find an increase in bias for the most 
luminous 10% of quasars. 

Some studies used photometrically detected galaxies 
rather than spectroscopic galaxies to increase the size of 
the tracer population and further dec rease their statisti¬ 
cal errors. Padmanabhan et al. (2009[) found no luminos¬ 
ity evolution in the cross-correlation between SDSS DR5 
quasars and luminous red galaxies with photom etric red- 
shifts at 0.2 < z < 0.6. jzhang et al. (2013) found no 
luminosity dependence in the cross-correlation between 
DR5 quasars and photometrically detected galaxies from 
SDSS Stripe 82. 

While most studies focused on the optical luminos¬ 
ity dependence of quasar clustering, other studies mea¬ 
sured th e luminosity depend ence at different wave¬ 


lengths. Hickox et al. (2009) found significant varia¬ 


tion in clustering strength between radio. X-ray, and in¬ 
frared se lected AGNs. Using X -ray selected quasars from 
ROSAT, Krumpe et al. (2012) found a 2 cr dependence of 
clustering strength on X-ray luminosity, and attributed 
the discrepancy between measurements of the optical lu¬ 
minosity dependence and measurements of the X-ray lu¬ 
minosity dependence to the much larger dynamic range 
in X-ray luminosity. 

Fewer studies have attempted to measure black hole 
mass dependent clustering due to the much higher uncer¬ 
tainties of black hole mass estimates. Black hole masses 
can be estimated from single-epoch spectroscop y using 
continuum luminosity an d emission line width (Vester- 
gaard & Peterson] 2006) with uncertaint ies up to 0.5 
decades for individual objects ( Shen][2M3 ). Throughout 
this paper, we refer to these estimates as “virial black hole 
masses” to emphasize the lar ge associated uncerta i nties. 

Using the same sample as Croom et al. (2005), Fine 
et al. (2006) constructed composite spectra tor each of 
10 bins in redshift and measured the virial mass for each 
composite spectrum. They found a Mbh ~ M dmh rela¬ 
tionsh ip in good agreement with the models of Ferrarese 
(2002), although with substantial uncertainties in the 
slope and zeropoint of the relation due to both the small 
dynamic range in virial mass and the Malmquist bias 
arising from the 2QZ flux limit. Other studies measured 
the e volution of c l usteri ng with virial mass at fixed red¬ 
shift: Shen et al. (2009) found no dependence of quasar 


clustering on virial mass using two bins in viral mass, and 
Zhang et al. (2013) found a 1-2 a difference in clustering 
strength between two bins in virial mass. However, these 
studies were hampered by a small number of bins and a 
small dynamic range in virial mass, implying an even 
smaller range in true mass due t o the uncertainty in the 


virial masses (Shen et al. 2009]). Komiya et al.] (2013) 


with larger samples and smaller error bars (da Angela 


studied the virial mass dependence ot AOIN clustering 
across a wide range of redshifts (0.1 < z < 1.0) and lu¬ 
minosities, combining samples from SDSS DR4 and DR7 
to measure clustering across 2 decades in virial mass. 
They found a significant trend of increasing clustering 
strength with increasing virial hole mass across 4 bins in 
virial mass, with a 2-3 cr difference in clustering strength 
between the highest and lowest mass groups. 

In this study, we measure the angular overdensity of 
photometric galaxies about spectroscopic quasars to ob¬ 
tain as numerous a tracer population as possible, min- 
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imizing statistical errors. We measure the luminosity 
and black hole mass dependence of the quasar-galaxy 
clustering amplitude at 0.65 < z < 0.9, using a galaxy 
sample drawn from the W ide-Field Infrared Survey Ex¬ 
plorer (Wrig ht et al.||2010| and quasar samples from the 
SDSS UR7 (Ishen et al. 1 1201 IB and DRIO quasar cata¬ 


logs (Paris et al. 


2014). Recause the WISE galaxies lack 


spectroscopic redshitts, we measure the angular correla¬ 
tion using only the galaxy and quasar positions on the 
sky. Since emission from galaxies at z ^ 0.8 peaks in the 
near-infrared, we expect that WISE selection will maxi¬ 
mize the size of our tracer population. 

In this paper we begin by discussing the quasar and 
galaxy samples and their selection criteria (Section . 
Then we present the angular clustering measurement 
(Section [3) and find the dependence upon luminosity and 
virial black hole mass (Section®. Finally we discuss the 
significance of these results ana compare them to previ¬ 
ous studies (Section]^. Throughout this paper, we use 
a A CDM cosmology wit h h = 0.7, Da = 0.7 and 17^ = 
0.3 (Spergel et al.| 200^ , matching the cosmology used 


to construct the DR7 and DRIO quasar catalogs. When 
computing the power spectrum of the linear-regime mat¬ 
ter field, we use the transfer function of |Eisenstein &:| 
|Hul ([19^ and U,, = 0.044, = 2.726 K, = 3, 

and ris = 0.93 (Spergel et al. 2003). All distances are 
measured in comoving h~^ Mpc. 

2. QUASAR AND GALAXY SAMPLE SELECTION 

We use a color-selected sample of z > 0.6 galaxies 
from WISE and 0.65 < z < 0.9 quasars from the SDSS 
DR7 and DRIO quasar catalogs to measure quasar-galaxy 
clustering. The galaxy sample was chosen to maximize 
purity rather than completeness, allowing us to obtain a 
high signal-to-noise measurement of quasar-galaxy clus¬ 
tering. Galaxy colors are measured by comparing SDSS 
and WISE imaging. 

2.1. Galaxy selection: SDSS and WISE 

The Wide-Field Infrared Survey Explorer is a satel¬ 
lite that mapped the entire sky at 3.4, 4.6, 12 and 22 
/rm at sensitivities corresponding to Vega magnitudes 
of 16.5, 15.5, 11.2, and 7.9 in unconfused regions of 
the sky (Wright et al. 2010). WISE collected at least 
12 exposures at each point on the sky, with coverage 
depth increasing rapidly towards the Ecliptic poles due 
to the WISE scan strategy. We select objects from 
the All-Sky Data Release Source Catalog, which con¬ 
tains 563 million objects that are not flagged as an im¬ 
age artifact and have both SNR > 5 in at least one 
band and detections in at least 5 single-band exposures. 
The All-Sky Data Release uses imaging taken from Jan¬ 
uary to August 2010, and extensive documentation can 
be found at http;//wise2.ipac.caltech.edu/docs/ 
release/allsky/. We use imaging from the W1 band 
at 3.4 /xm, which has angular resolution of 6" (Wright 
et al.|20l0 1. The WISE limiting magnitude varies across 
the sky, most notably due to increased source confusion 
at low Galactic latitud^ We only consider objects with 
Galactic latitude b > 25°; at these latitudes, the W1 

^ See Figure 9 at http;//wise2.ipac.caltech.edu/docs/ 
release/allsky/expsup/sec2_2.html 


magnitude limit is approximately constant. WISE pho¬ 
tometry is given in Vega magnitudes, uncorrected for 
Galactic dust extinction, which is negligible at 3.4 /tm. 
For sources with 15.5 < W1 < 16, similar in brightness 
to those selected in our galaxy catalog, the astrometric 
accuracy of WISE is « 0.4", as measured by the RMS of 
the WISE-2MASS positional differenc^ 

The Sloan Digital Sky Survey imagecUan ar ea of 14555 
deg^ in 5 filters {ugriz-, Fukugita et al. ]|1996|) mostly at 
high Galactic latitude, using a 2.5 m wide-tie ld telescope 


at Ap ache Point Observatory in New Mexico (Gun n et al. 
200611 and a camera with 30 2048 by 2048 GGDs ( |Gunn 


et al.|1998 1. Imaging data was collected only under pho- 


tometric conditions (Hogg et al. 2001) and photometry 

was calibrated to 1% accuracy 

(ISmith et al.||2002 

Ivezic 

et al.|2004 

Tucker et al.|2006 

Padmanabhan et al. 

2008| 


with negligible spatial vari ation in photometric calibra- 
tion (Fukugita et al. 2004|). All SDSS magnitudes cited 
in this paper are corrected for Galactic e xtinction using 


the S chlegel-Finkbeiner-Davis dust map (Schleg el et al. 
1998t and are reported as asinh magnitudes ([Luptm 


et al. 1999) in the AB system (Fukugita et al. 1R9^~ 


SPSS ^oto metry is 95% complete at r = 22 (Stoughton 
et al.||2002|). Astrometry is typically accurate to 0.1" 


(|Pier et al.|2003|, although the photom etric data used 
in this paper (DR8; Aihara et al. 2011b I contain an as¬ 


trometric calibration error that causes a shift of 0.24" 
north and 0.05" west over a large re gion covering most o f 
the survey with declination > 41° ( Aihara et al.|[Mlla l. 
However, the effect of the astrometry error on our work 
is negligible. 

The SDSS photometric pipeline separates extended 
sources from point sources based on the difference be¬ 
tween the PSF magnitude and a composite model mag¬ 
nitude consisting of a linear combina tion of de Vau- 
coule urs and exponential light profiles (Stoughton et al. 
2002). Many galaxies at z ~ 0.8 are unresolved, par¬ 
ticularly in regions of poor seeing. To eliminate seeing- 
dependent variations in galaxy sample density, our sam¬ 
ple contains point sources as well as extended sources 
in SDSS imaging. For both point sources and extended 
sources, we use the best-fit exponential or de Vaucouleurs 
model magnitudes (Stoughton et al.|2002 1. We select our 
objects using the catalogs from the DR8 data sweeps, 
which contain point sources with at least one extinction- 
corrected PSF magnitude less than [22.5, 22.5, 22.5, 22, 
21.5] {ugriz) and extended sources with at least one 
extinction-corrected m odel magnitude less than [21, 22, 
22, 20.5, 20.1[ {ugriz) ( Blanton et al.|[^005 1. 

We use a color cut to ensure that our sample is com¬ 
posed primarily of z > 0.6 galaxies. In order to choose 
an appropriate color cut, we match WISE and SDSS 
photometry to spectroscopic redshifts for 20,000 galax¬ 
ies measu red by the AGN and G alaxy Evolution Survey 
(AGES) (Kochanek et al. 2012). To calculate WISE- 
SDSS colors, we use r magiiitudes in the AB system and 
W1 magnitudes in the Vega system. Golor cuts based 
on WISE bands alone are inadequate, so instead we se¬ 
lect all galaxies with r — W1 > 5.5. This is a conser¬ 
vative color cut that minimizes the number of stars and 
low-redshift galaxies in our sample, reducing statistical 

^ http://wise2.ipac.caltech.edu/docs/release/allsky/ 
expsup/sec6_4.html 
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uncertainty in our clustering measurement. As Figure 
shows, the vast majority of the galaxies satisfying our 
color cut have z > 0.65. However, many galaxies with 
z > 0.65 have r — W1 < 5.5, particularly blue galaxies 
with 0.6 < z < 0.8. In this study, we assume that the lu¬ 
minosity and virial mass dependence of clustering is the 
same for red galaxies as for blue galaxies. 


8 -1-1-1-1-1—I—I-1-1-1-1-1-1—I—r 



Q _I_I_I_I_I_1-1-1-1_I_I_I_I_I-1- 

0 0.2 0.4 0.6 0.8 

Redshift 

Figure 1. Redshift vs. extinction-corrected color for a sample 
of galaxies from the AGES survey. Note that r is reported 
in AB magnitudes while lEl is reported in Vega magnitudes. 
All points above the red line would be included by our color 
cut. Since most of these galaxies have z > 0.65, our color cut 
removes many contaminating z < 0.6 galaxies. Almost all of 
the severe deviations from the linear trend are objects with 
spurious colors and should be ignored. 

We begin by choosing WISE detections with 
14 < WI < 16.5, W1 SNR > 5 and Galactic latitude 
b > 25°. Less than 0.1% of the BOSS luminous red 
galaxies with 0.65 < z < 0.9 have W1 < 1'^ so cut¬ 
ting objects with W1 < 14 reduces stellar contamina¬ 
tion without removing z > 0.6 galaxies. We also remove 
potentially variable, saturated, or contaminated detec¬ 
tions and detections that were deblended more than once 
(WISAT > 0, W1CC_MAP > 0, VAR_FLG > 5, and 
NB > Since nearly every z > 0.6 galaxy is unre¬ 
solved in \VISE imaging, we remove detections that are 
extended sources in WISE (EXT_FLG > 0). Last, we 
remove detections with more than 10% of observations 
contaminated by scattered moonlight (MOON_LEV > 
1). Since the detections with moonlight contamination 
are highly clustered due to the WISE moon avoidance 
maneuver^ we exclude all WISE detections in regions 
with high moon contamination. 

Visual inspection of WISE detections reveals a prob¬ 
lem with WISE deblending. About 5-10% of WISE de- 

^ These are mostly CMASS galaxies; see|White et al.|(|201lf for 
the sample selection criteria. 

^ See http://wise2.ipac.caltech.edu/docs/release/allsky/ 
expsup/sec2_2a.html for definitions. 

See http://wise2.ipac.caltech.edu/docs/release/prelim/ 
expsup/sec3_4a.html 


tections within the SDSS imaging footprint but with¬ 
out a matching SDSS detection are blended into bright 
neighboring detections in WISE imaging (see Figure 0. 
In many of these cases, SDSS found a bright extended 
galaxy within 10" of the WISE detection. It appears 
that the WISE detection should have been matched to 
the SDSS galaxy, but differences in the WISE and SDSS 
deblending algorithms caused a spurious offset between 
the WISE and SDSS positions, large enough so that the 
two detections could not be matched. The SDSS galaxies 
are all contaminating z < 0.6 galaxies that should not be 
included in our sample. Note that we still find many of 
these deblending problems even when we only consider 
WISE detections with NB < 2. 



8" aperture minus fitted magnitude 

Figure 2. Top: example of a WISE blending prob¬ 
lem in WISE (left) and SDSS (right), located at 
22’'52™9°.02,1°38'6.4". The two images cover the same re¬ 
gion of the sky. The colored circles on the WISE image are 
centered on WISE detections; the red circle lacks an SDSS 
match. It is clear that the red and blue circles are both asso¬ 
ciated with the bright yellow galaxy at the center of the SDSS 
image and the red circle is merely an artifact of the WISE de¬ 
blending process. Bottom: Distribution of Wl aperture mi¬ 
nus fitted magnitude (A — F) for WISE detections included 
in our sample. The tail on the left results from deblended 
sources such as those in the top images, which have an ex¬ 
tended flux profile despite being classified as point sources. 
Because these could be contaminating low-redshift objects, 
we removed all detections with A — F < —0.5. 

In order to remove these detections from our sample, 
we compare two different measurements of W1 magni¬ 
tude: one computed by summing all the flux within 
an 8.25" aperture {W1MAG_2] referred to as A for 
clarity), and the standard Gaussian profile-fitted mag¬ 
nitude measurement {WIMPRO; F for clarity). We 
expect detections with blending problems to have unusu¬ 
ally extended flux profiles. Since the profile-fitted mag- 
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nitudes are largely determined by the central flux, detec¬ 
tions with blending problems should have brighter 8.25” 
aperture magnitudes than proflle-fltted magnitudes. In¬ 
deed, the distribution of aperture minus fitted magni¬ 
tudes {A — F) has an asymmetrical tail of detections with 
an extended flux profile (Figure]^. Since the distribution 
is nearly zero for A — F > 0.5, we impose a symmetrical 
cut on the opposite tail and keep only those detections 
with A — F > —0.5 (Figure]^. Visual inspection of 50 
objects with A — F < —0.5 shows that nearly three quar¬ 
ters have blending issues like Figure 1^ However, blend¬ 
ing issues are very uncommon for A~— F > —0.5. We 
And that cutting all objects with A — F < —0.5 is more 
effective at removing deblending problems than cutting 
all detections that were deblended at least once. 

Our galaxy sample consists of all WISE detections 
meeting the above criteria that have r — W1 > 5.5. 
Since SDSS photometry is complete to r = 22 and our 
sample only includes detections with W1 < 16.5, WISE 
detections imaged by SDSS but lacking an SDSS match 
meet our color cut and are included in our sample. Our 
galaxy sample consists of two components: WISE-SDSS 
matches with r — W1 > 5.5 (“matches”) and WISE de¬ 
tections located within the SDSS imaging area that do 
not match an SDSS detection (“non-matches”). In order 
to determine the SDSS imaging area, we use the SDSS 
imaging mask of Ho et al. (2012|, which is more restric¬ 
tive than the SDSS imagihg footprint: it removes re¬ 
gions with poor seeing, low Galactic latitude, and nearby 
bright stars, all of which may contaminate clustering 
measurements. This yields a Anal area of 6966 deg^. 
We exclude all W ISE detect i ons ly ing outside the SDSS 
imaging mask of Ho et al. (20121 because we have no 
information about these objects’ VVISE-SDSS colors. 

We apply additional cuts to the WISE-SDSS matches 
to ensure that the SDSS detections are real objects rather 
than artifacts. We remove duplicate d etections by only 
includ ing “survey primary” detections (Stoughton et al. 
20021. For matches with r < 22, we remove all de¬ 


tections with SDSS imaging flags indicating dirty pho¬ 
tometry (e.g. saturated detections, cosmic ray strikes, 
image-processing artifacts, etc|^l. We do not remove de¬ 
tections with r > 22 and dirty ^otometry because these 
detections are faint enough to be spuriously flagged as 
“dirty” even if they are real objects. We also remove all 
detections with r-band dust extinction greater than 0.3 
magnitudes. 

Ignoring the WISE resolution of 6", the astrometric 
precision of SDSS and WISE implies that all WISE-SDSS 
matches should have a separation < 1". However, since 
WISE has a much larger angular resolution than SDSS, 
we And many matches with separations between 1" and 
3". Visual inspection of SDSS and WISE images shows 
that many matches with separations > 1" are not merely 
a result of imprecise WISE astrometry. Instead, two 
SDSS detections separated by < 6" are merged into a 
single WISE detection in a substantial fraction of these 
matches. If the second-closest SDSS detection is sub¬ 
stantially brighter than the match, it is possible that a 
large portion of flux from the WISE detection should 
have been assigned to the second-closest SDSS detection 

® See https://www.sdss3.org/dr8/algorithms/photo_flags_ 
recommend.php tor detinition. 


rather than to the match. As a result, the match may not 
actually meet our color cut. We remove these problem¬ 
atic matches by defining the SDSS match as the brightest 
SDSS detection within 3" of the WISE detection. This 
definition removes matches in which the second-closest 
SDSS detection is brighter than the match, leading to 
a spuriously high r — W1 for that WISE detection. By 
excluding these detections, we reduce contamination in 
our sample. 

Our Anal galaxy sample contains 4,168,855 objects, in¬ 
cluding 2,675,189 non-matches and 1,493,666 matches 
(composed of 331,327 point sources and 1,162,339 ex¬ 
tended sources in SDSS imaging). We estimate the sever¬ 
ity of contamination from Galactic stars by plotting the 
sky distribution of our sample in Galactic coordinates. 
We also plot the sky distribution of our sample in Ecliptic 
coordinates to determine whether our sample is affected 
by fluctuations in the WISE coverage depth due to over¬ 
lapping scan^(Figure]^. We observe substantial gradi¬ 
ents on the sky related to the Galaxy but do not observe 
significant gradients following the WISE scan pattern. 
Both the SDSS-identifled extended sources and the non¬ 
matches have densities ~80% greater at high Galactic 
latitudes. The anti-correlation with the Galactic center 
results from source confusion: at low Galactic latitudes, 
a galaxy is more likely to be masked by a Galactic star. 
The similarity between the sky distributions of the ex¬ 
tended sources and the non-matches suggests that the 
non-matches are largely composed of z > 0.6 galaxies 
rather than Galactic stars or spurious WISE detections. 

SDSS-identifled point sources display the opposite pat¬ 
tern, with densities nearly twice as high at low Galactic 
latitudes than at the poles. This indicates that a sub¬ 
stantial fraction of the point sources are very red Galac¬ 
tic stars. Nevertheless, we believe it is prudent to retain 
the point sources in our sample. The SDSS star/galaxy 
separation is seeing dependent, with galaxy density de¬ 
creasing in ima ging with poor seeing , particularly for 
faint detections (Scranton et al. 20021. As a result, ex¬ 
cluding the point sources from our sample would lead to 
seeing-dependent spatial variations in density, which are 
likely correlated with seeing-dependent spatial variations 
in quasar density. Given that point sources compose only 
7% of our sample and that the point sources correlate sig¬ 
nificantly with the quasars (see Figure]^, we believe that 
is appropriate to keep the point sources in our sample to 
eliminate the possibility of seeing-dependent systematic 
error. 


2.2. Quasar selection 

In order to obtain as large a luminosity range as possi¬ 
ble, we measure the quasar-galaxy clustering amplitude 
for quasars selected from two quas ar catalogs, the SDSS 
DR7 catalo g with 105,783 qu asars (Shen et al. 2011) from 
SDSS I/II (lYork et al.|2^, and the DRIP catal og with 
166,583 quasars (Paris et al. 2014|) fro m BOSS (Eisen- 
stein et al. 2011 Dawson et al. 2013|. In both SDSS 
1/11 and BUSS, quasar candidates were selected from 
SDSS photometry using object colors. SDSS I/H tar - 
geted objects with 15 < i < 19.1 (Schneider et al.||2010|. 


^ The WISE scan pattern in is given in Figure 5 in 
http://wise2.ipac.caltech.edu/docs/release/allsky/expsup/ 
sec6_2.html 
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Non-matches 



SDSS extended sources 



SDSS extended sources 




SDSS point sources 



SDSS point sources 



Galaxy Density 


466 _ 


7t2 deg' 


DR7 Quasars 


DR10 Quasars 



Figure 3. Top and middle (left to right): Densities of the entire galrixy sample; all WISE galaxies lacking an SDSS match; 
all WISE galaxies matched to an SDSS extended source; and all WISE galaxies matched t o an SDSS point source. Both plots 
only display HEALPix pixels with more than 20% of their area covered by the|Ho et al. (20121 DR8 imaging footprint. Top: 
Ecliptic coordinates with 180° Ecliptic longitude running through the center or the plot and longitude increasing to the left. 
Middle: Galactic coordinates with Galactic longitude increasing clockwise from the Galactic center at the bottom of each plot. 
While the variation in density with Galactic latitude is clearly apparent, we do not observe any variation in density in regions 
of deeper WISE coverage (wide strip between 200° and 230° Ecliptic longitude and -90° to 90° Ecliptic latitude). Bottom: 
Comparison of galaxy and quasar sky distribution. Bottom left: galaxy density used in the calculation of the quasar-galaxy 
angular overdensity. Pixels with less than 20% of their area covered by the imaging footprint are excluded. Bottom center and 
right: DR7 and DRIO quasar distribution. Because of the low density of quasars, the fluctuations between pixels are dominated 
by shot noise rather than by intrinsic variation in quasar density. Only quasars in the region of sky covered by both DR7 and 
DRIO spectroscopy are displayed. 
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while BOSS targeted ob jects with eit her r < 21.85 o r 
g <22 (Paris et al. 2014) and i > 17.8 (Bovy et al.|2011 1. 
Quasar candidates w ere arranged on spectroscopic plates 
(Blanton et al. ]200^ , o bserved using a twi n multi-object 
iiber-ted spectrograph (Smee et al. 2013), and reduced 


mates assume that the quasar’s broad-line region (BLR) 
is virialized: 


Menc = A7 bH = 


^BLR-^BLR 

G 


( 2 ) 


and classified as a quasar, galaxy or star (Bolton et al. 


20121. The wavelength range is 3800 to 9200 A for SD^ 


I/II and 3600 to 10400 A for BOSS ( Smee et al.||2013 1. 


DR7 quasar candidates were selected using a uniform 


The BLR velocity is inferred from the width of a par¬ 
ticular broad line. The radius is determined by measur¬ 
ing the continuum luminosity: reverberation mapping of 
z < 0.3 AGN found a tight relatio nship between con - 


targe t algorithm based on their colors (Richards et al. tinuum luminosity a,nd BLR radius ( Kaspi et al. 


2000 ). 


imaged objects selected by other targeting algorithms, 
and we reject the half of the quasar catalog that was not 
selected using the uniform targeting algorithm. Quasars 
in DRIO are selected using a variety of methods. Half 
are selected uniformly to form the COR E sample, using 
the XDQSO method (Bovy et al. 2011), and the other 
half are selected inhomogeneously to maximize surface 


density (see Paris et al. 
selection algorithms). 


2014 for details on the various 
\Ve only use the CORE sample 


in this work. Both the DR7 uniform quasars and the 
DRIO CORE quasars are uniformly distributed across 
the sky. The DR7 spectroscopic footprint is substantially 
larger than the DRIO footprint, because DR7 is the final 
data release for SDSS II while DRIO is not the final data 
release for SDSS III. We only use quasars that lie in the 
intersection of the DR7 and DRIO footprints. We also 
only consid er quasars lying within the DR8 photometric 
footprint of Ho et al. (2012), since our galaxy sample is 
restricted to this footprint. 

The final quasar sample contains 7,049 quasars with 
0.65 < z < 0.9, 4,206 from DRIO and 2,843 from DR7. 
Figure shows the redshift, luminosity and virial mass 
distributions for the DR7 and DRIO quasars. To measure 
the luminosity dependence of the clustering amplitude, 
we split the DR7 quasars into three groups by luminosity 
and the DRIO quasars into four groups by luminosity. 
Similarly, we split the DR7 quasars into four groups by 
virial mass and the DRIO quasars into four groups by 
virial mass. 

Redshifts in both t he DR7 and DRIO catalogs are ac- 


2002). The DR7 quasar catalog includes serendipitously This relationship arises because the size of the BLR is 

regulated by the amount of ionizing radiation emitted by 
the quasar, which is proportional to the optical contin¬ 
uum luminosity. Since the reverberation mapping sam¬ 
ples used to calibrate the single-epoch virial mass esti¬ 
mates use Ha and H/3 line widths, viri al mass esti mators 
based on Ha and H/3 are most reliable (Shen 2013). How¬ 
ever, H/3 is redshifted beyond the edge of the DR7 spec¬ 
trograph (9200 A) at z > 0.85, so we instead use Mgll 
based mass estimates . Mgll line widths correlate well 
with H,5 line widths (Shen 2013), but the Mgll masses 
may possess substantial systematic errors due to the lack 
of Mgll-based reverberation mapping masses. However, 
since we are interested in the slope rather than the nor¬ 
malization of the black hole mass-clustering strength re¬ 
lationship, we are not concerned with systematic errors 
resulting in a uniform offset in black hole mass. 

We are primarily concerned with two kinds of error 
in the virial mass estimates: increased scatter due to 
the large uncertainties in the virial mass estimates and a 
luminosity-dependent bias arising from the flux limit of 
the sample. First, the spread in virial mass is wider than 
the spread in true mass because of the large uncert ain- 
ties in individual virial mass estimates. Shen] (2013) es¬ 
timated that the scatter in virial mass at fixed true mass 
is up to « 0.5 decades, arising from both measurement 
errors in line width and continuum luminosity, and un¬ 
certainties in the virial mass calibrations. We reduce the 
measurement uncertainties slightly by excluding quasars 
with poor continuum or Mgll emission line fits (y^/d.o.f 
> 2), but this does not reduce the calibration uncertain¬ 
ties, which are the dominant source of scatter in virial 
mass at fixed true mass. Because of the large uncer¬ 
tainty in virial mass, the dynamic range of true mass in 
our sample is less than the measured range in virial mass, 
and the mean true mass for quasars in a given virial mass 
bin is less extreme than the mean virial mass. 


curate to Az < 0.01 (Schneider et al. 2010, Paris et al. 


2014). In this paper, we measure luiinnosity using the 
absolute i-band magnitude K-correct ed to z = 2 for both 
DR7 and DRIO (Richards et al.|2006 ). These magnitudes 
measure the quasar luminosity in a bandpass centered at 
2500 A in the rest frame. We convert absolute magni¬ 
tude to 2500 A luminosity in erg s“^ using Equation 4 
from Richards et al. (2006) and then to solar luminosity 
using Lq = 3.827 x 10'’’'’* erg s“^. Reported i-band mag¬ 
nitude errors for quasars in our sample are < 0.03 for 
DR7 and <0.1 for DRIO. Note that the DR7 quasars 
are substantially more luminous than the DRIO quasars 
(Figure]^, since BOSS targeted fainter, higher-redshift 
quasars compared to SDSS I/H. 

We use the sin gle-epoch virial mass estimates from 


Shen (2013) also discussed a luminosity-dependent bias 


Shen et al. (2011|) for DR7 quasars, and estimates com- 
puted using similar methodology for the DRIO quasars 
(Y. Shen, private communication) The virial mass esti- 


^ The line fitting methodoloKy for DRIO is described in |Shen| 
I&; Liu| (I 2 OI 2 I 1 and at http: //users. obs. carnegiescience. edu/ 
yshen/Bh_mass/dr9 .htm i'he practical impact of the ditterent lit- 
ting procedures tor DK7 and DRIO is negligible (Y. Shen, private 


arising from uncertainties in line width and luminosity. 
Measurement errors, scatter in the radius-luminosity re¬ 
lationship, non-virial motion, and a time lag between 
changes in luminosity and radius may lead to uncorre¬ 
lated errors in luminosity and line width. If uncorrelated 
errors are present, virial masses in flux-limited samples 
will be biased high relative to true masses, since the in¬ 
crease in average luminosity caused by the flux limit will 
not be entirely cancelled by a decrease in line width. The 
magnitude of the bias will be greatest for low luminosity 
subsamples that lose a subst antial number of qua s ars to 
the flux limit. Shen (2013) and Shen & Kelly (2012[ 
found evidence tor luminosity-dependent bias in IVlgl 
masses. Since less massive quasars are also less lumi- 

communication). 
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nous, the luminosity-dependent bias will be stronger for 
low mass quasars than for high mass quasars. However, 
since each bin in virial mass covers a large range in true 
mass, the magnitude of the luminosity-dependent bias is 
quite similar for all virial mass bins. Because our quasar 
samples are flux limited, the luminosity-dependent bias 
causes the virial masses to systematically overestimate 
the true masses, but this uniform offset in black hole 
mass does not affect our measurement of the slope of the 
black hole mass-clustering strength relationship. 

2.3. Systematic effects 

The difference in redshift distribution between different 
quasar groups leads to systematic differences in quasar- 
galaxy clustering strength. The quasar l uminosity func¬ 
tion evolves significantly with redshift (Richards et al. 
20061, so high luminosity quasars have a higher mean 
redshift than low luminosity quasars. Angular clustering 
strength varies with redshift due to both r edshift evolu¬ 
tion o f the quasar autocorrelation function (Groom et al. 
20051 and the redshift distribution of our galaxy sample. 
Therefore, any difference in clustering between quasar 
groups with different luminosities may in fact arise from 
the difference in redshifts. 

To isolate the luminosity dependence of clustering, we 
assign a redshift-dependent weight to each quasar to force 
the weighted redshift distributions of each subsample to 
match the DRIO redshift distribution. For each subsam¬ 
ple, we place the quasars into bins of width Az = 0.01 
and compute the weights by dividing the DRIO redshift 
distribution by the subsample’s redshift distribution. To 
ensure that no quasars are weighted by more than 3 or 
less than 1/3, we combine the two lowest-luminosity DR7 
groups (see Figured. 

Because the densrfy of both the WISE galaxies and the 
stellar contaminants in the WISE sample vary across the 
sky, differences in large-scale sky distribution may lead to 
differences in clustering amplitude. Using only quasars 
located within the DR7/DRIO overlap region largely alle¬ 
viates this problem by forcing the sky distribution of DR7 
and DRIO quasars to match (see Figure 1^. While the 
resulting sky distributions are not identical, the residual 
variation in both WISE galaxy density and contaminat¬ 
ing star density is quite small. We eliminate the effects 
of residual variation in WISE galaxy density by measur¬ 
ing t he galaxy density sep arately in each HEALPix pixel 
(see Gorski et al. (20051 for details about HEALPix). 
The stellar contaiinnation, as measured by the average 
fraction of SDSS point sources in the HEALPix pixel 
surrounding each quasar, is just 1.2% greater for DRIO 
quasars than for DR7 quasars. The impact of varying 
stellar contamination is therefore considerably lower than 
the magnitude of the measured overdensity w. 

Small scale variations in galaxy density may also af¬ 
fect our clustering measurements. While Figure can¬ 
not display density variations on scales less than a few 
degrees, we expect to observe variations in WISE den¬ 
sity on arcminute scales due to both source suppression 
near bright sources and image artifacts within the ha¬ 
los of very bright source^ The vast majority of these 
bright sources are uncorrelated with quasar positions, 

® See http://wise2.ipac.caltech.edu/docs/release/allsky/ 
expsup/sec6_2.html 


SO we do not apply a mask to remove the area around 
bright WISE sources. Unlike the large-scale variations 
in quasar and galaxy densities, small-scale variations in 
galaxy density should have the same effect on both DR7 
and DRIO groups, and on different groups in luminos¬ 
ity and virial mass. However, bright quasars in WISE 
imaging may lead to both an overdensity of nearby im¬ 
age artifacts and an underdensity of nearby faint sources 
due to the increased background. While we failed to And 
image artifacts around the 50 brightest quasars in Wl, 
these quasars, with Wl « 12, suppress the density of 
Wl = 16.5 sources at separations < 18 '(3 As a result, 
we restrict our measurement of the quasar-galaxy cross¬ 
correlation function to angular scales > 18", correspond¬ 
ing to « 0.2 h~^ Mpc. At these scales, we expect no sig¬ 
nificant suppression of galaxy density due to increased 
background from nearby quasars. Similar source sup¬ 
pression in SDSS is only observed at separations less than 
15" for galaxies with similar brightness as our galaxy 
sample about stars with s imilar brightness as our quasar 
sample (Ross et al.||2011 1. 


3. MEASURING THE ANGULAR OVERDENSITY 
3.1. Methods 

In this paper, we measure the angular overdensity 
w{rp) rather than the three-dimensional correlation func¬ 
tion ^(r) or the projected cross-correlation function 
Wp{rp). The three-dimensional quasar-galaxy c ross¬ 
correlation fun ction is well-fit by a power law (e.g. Yee 
fc Green|[l9^ : 

«*■) = ( 7 )’ (3) 

where rg is the correlation length, which is typically used 
to characterize clustering strength. The projected cross¬ 
correlation function Wp{rp) is the integral of the three- 
dimensional cross-correlation function along the line of 
sight: 


/ OO 

^(rp, tt) dTT = 

-00 


r(i/2)r[(7-i)/2]^i_^ 

P 


r(7/2) 


.(4) 

where T is the gam ma function and tt is the line-of-sight 
distance. Following Zhang et al. (20131, we can relate the 
projected cross-correlation function Wp(rp) to the angu¬ 
lar overdensity w{rp): 


Wpirp) = { — / 


(5) 


where n is the WISE density near each quasar (Figure 
and po is the density of galaxies at each quasar’s redshut 
that meet our color and Wl cuts. 

Estimating rg and bg requires finding the projected 
cross-correlation function, which requires an estimate of 
pg. One can compute pg by integrati ng the galaxy lumi¬ 


nosity function at ea ch redshift (e.g. Zhang et al. 2013 
Komiya et al. 20131. However, this method is insuffi¬ 


ciently precise to measure the host halo mass to a reason¬ 
able accuracy: since bias is a shallow function of host halo 
mass at z « 0.8, small errors in the luminosity function 

See Figure 25 at http://wise2.ipac.caltech.edu/docs/ 
release/allsky/expsup/sec6_2.html 
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Figure 4. Left: redshift distribution for DR7 and DRIO quasars. Center-. Distributions of Mi{z = 2) for DR7 and DRIO 
quasars (« 2500 A restframe wavelength) in solar luminosities. The vertical lines divide the sample Into 4 DRIO groups and 3 
DR7 groups. The right-most cut in the DRIO distribution at 11.43 log L/Lq minimizes the overlap between the most luminous 
DRIO sample and the least luminous DR7 sample. Right-. Virial black hole masses for DR7 and DRIO quasars, computed by 


Shen et al. (20111 from continuum luminosity and emission line widths. The DR7 and DRIO quasars are each divided into 4 
groups; the divisions are not shown because several of the bins overlap. 


(and thus pq, vq, and the bias) lead to relatively large er¬ 
rors in the host halo mass. For instance, small differences 
in photometric calibration between SDSS and the instru¬ 
ment used to find the luminosity function will lead to sub¬ 
stantial errors in the luminosity function, causing large 
uncertainties in the host halo mass. Moreover, to find po 
we would need to apply our color cut to the galaxy lu¬ 
minosity functions, requiring knowledge of galaxy SEDs 
at z « 0.8. Ultimately, finding tq and 5 q is not neces¬ 
sary to measure the change in clustering with luminosity 
and virial mass, since such changes will be just as appar¬ 
ent in the angular overdensity w(rp). Given the difficul¬ 
ties associated with computing tq and bq from our data, 
we choose to only measure the clustering dependence of 
w{rp) in this paper. 

We measure the angular quasar-galaxy overdensity by 
counting the number of excess galaxies at angular sep¬ 
aration 9 from each quasar. We count the number of 
galaxies N in an annulus of width A9 and divide by the 
expected density of galaxies: 

»(")= (E 5;)/(E ”'*)-! (6) 

where n is the average galaxy density, A is the area in¬ 
side the SDSS imaging footprint, the index i ranges over 
every quasar, and the weights w are computed from the 
redshift distribution as described in Section [231 Because 
the galaxy density varies substantially across the sky, n is 
estimated by using the galaxy density in each HEALPix 
pixel (Figure]^. We exclude 18 quasars lying in pixels 
with less than 20% coverage of the imaging footprint; 
for these pixels, Poisson variations in the density are > 
1%, comparable in magnitude to the Galactic gradients 
shown in Figure and thus we consider the density es¬ 
timates for these quasars to be unreliable. 

We determine A, the area within the imaging mask, 
using a Monte Garlo method. We create a catalog of 
random points lying outside the imaging footprint with 
density n-random 2909 deg“^ much larger than the maxi¬ 
mum n, 712 deg“^ (Figure^. We use this catalog to find 
the area within the imaging footprint. A, by subtracting 
the area outside the imaging footprint from the area of 


the full annulus, Afuip 

A = (7) 

^random 


where the area outside the imaging footprint is de¬ 
termined by dividing the number of random points in 
the annulus, N^andom-, by the density of random points 
Tirandom- Next we transform from w{9) to w{rp), where 
Tp, measured in comoving h~^ Mpc, is the distance 
between the quasar and the galaxy assuming that the 
galaxy lies at the same redshift as the quasar. We com¬ 
pute the area of the full annulus, Afuu, using Vp rather 
than 6: 


Afuii = tt 


t .2 _ r .2 


( 8 ) 


where ri = Tp — Avp, r 2 = Vp + Avp, and Da{z) is the 
angular-diameter distance to a quasar at redshift z. 

We measure w{rp) in the following 5 bins (in units of 
Mpc): 0.2-0.4, 0.4-0.8, 0.8-1.6, 1.6~3.2, and 3.2- 
6.4. We cannot measure w{rp) at separations less than 
0.2 h~^ Mpc (angular separation 18") because at smaller 
separations the density of lUl = 16.5 galaxies is sup¬ 
pressed by the wings of the central quasar’s flux profile. 
At angular scales larger than 6.4 h~^ Mpc, we find that 
our measurement o f w( r„) is contaminated by systematic 
errors (see Section 3.31. 

We use bootstrap resampling to estimate our error 
bars. We resample by HEALPix pixel: from the 154 pix¬ 
els with at least one quasar, we randomly select 154 pix¬ 
els with replacement and measure w(rp) for all quasars in 
the selected pixels. We use 50,000 resamples to calculate 
error bars and the covariance matrix. Equation gives 
the reduced covariance matrix for the measure ment of w 
across the entire sample, Rij = CijI^^/CuCjj, where C 
is the covariance matrix: 
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Since vjp{rp) is related by vj{rp) by a constant of pro- 
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portionality (Equation [^, we fit a power law to w{rp): 

w{rp) = Pr-^ ( 10 ) 

where /3 oc rj and d = 7 — 1. We use minimiza¬ 
tion to find the best fit values of /3 and 6 for our mea¬ 
surement of w{rp) using the entire quasar and galaxy 
samples. Since the off-diagonal terms of the covariance 
matrix are nonzero (Equation , we compute using 
the full covariance matrix. Using bootstrapping, we con¬ 
firm that the sampling distribution of w(r) in each bin 
is very well approximated by a Gaussian, indicating that 
X^ minimization is an appropriate curvefitting technique. 
We will use the clustering amplitude /3 to characterize 
the strength of quasar-galaxy clustering for each of our 
subsamples. 


3.2. Results 

Figure [^shows the measured angular overdensity and a 
best-fit curve, using a sample of 7,049 quasars from both 
DR7 and DRIO. A two parameter fit yields a minimum 
of 10.23 with 2 degrees of freedom, /3 = 0.105 ±0.007, 
and 6 — 0.84 ± 0.05. Varying the minimum bin radius 
does not substantially affect /3 or 6, nor does successively 
removing each bin from the fit. 

By dividing w{rp) by the projected correlation func¬ 
tion of the linear-regime matter field, we obtain a quan¬ 
tity proportional to 6 qq, the square of the linear quasar- 
galaxy bias (Equation [^. We plot this quantity in Fig¬ 
ure [ 5 ] Quasar-galaxy clustering at 0.2-6.4 h~^ Mpc 
arises from a mixture of one-halo and two-halo terms: 
the one-halo term refers to clustering within the same 
dark matter halo, while the two-halo term refers to clus¬ 
tering between different halos. One-halo clustering leads 
to an increa se in the linear bias at Vp < 1 h~^ Mpc (Shen 
et al.||2013) and cosmological hydrodynamic simulations 
indicate that the on e-halo term dominat es clustering at 
Up < 0.3 h~^ Mpc (Degraf et al. 20111. We interpret 
the sham increase in linear bias at Vp = 0.3 h~^ Mpc 
(Figure^ as evidence for one-halo clustering in this bin. 

3.3. Testing for systematics 

We check for systematic effects by measuring w{rp) for 
WISE galaxies about 37,402 quasars with 2.2 < z < 3.5. 
We randomly assign each quasar a redshift 0.65 < z < 
0.9 so that we measure w{rp) on the same angular scales 
as for our z ~ 0.8 quasar sample. We expect a small 
signal due to the gravitational lensing of high-redshift 
quasars by z ~ 0.8 galaxies. The galaxies magnify the 
high-redshift quasars, lowering the flux limit in the region 
near the ga l axy and creating a cross-cor r elation (Myers 
et al.||20^ |2005[ [Scranton et al.|[2005[ jlVlenard et al. 
2010||. The strength of this signal is proportional to the 
magnification /r: 


t \ ^ 

w\rp) oc fi ^ 


( 11 ) 


where 9e is the Einstein radius for a galaxy in our sam¬ 
ple (typically « 1” for galaxies at z « 0.8) and 9 is the 
angular separation between the galaxy and the quasar. 
The strength of this signal depends on both the level 
of contamination by low-redshift galaxies and stars and 
the average apparent magnitude of the quasar popula¬ 
tion. A linear least squares fit of constant Tpwfrp) yields 


Tpwfrp) = —0.0018 ± 0.0023. However, for bins centered 
at 1000" and 2000" (9.6 and 19.2 h~^ Mpc at z = 0.8, 
respectively), we find a 2-3 a deviation from zero, with 
w{rp) « 0.003 for these two bins. Since the gravitational 
lensing signal peaks at much smaller angular scales, we 
attribute this large-scale overdensity to a systematic er¬ 
ror present on all scales. While this error is very small 
compared to w{rp) at smaller scales, it is close enough to 
w{rp) for the bins at Vp = 9.6 and 19.2 Mpc that we 
restrict our measurement of angular clustering to sepa¬ 
rations of 0.2-6.4 h~^ Mpc. 
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Figure 6. Correlations between the DR7 and DRIO quasars 
and SDSS-identified extended sources (red squares, dotted 
line), non-matches (blue triangles, long-dashed line), and 
point sources (black circles, solid line). The lines are linear 
least squares curvefits using 5 — 0.84, the slope from the two 
parameter fit in Figure Angular clustering amplitudes at 1 
h~^ Mpc: extended sources, /3 = 0.134 ± 0.009; non-matches, 
P = 0.098 ± 0.008; point sources, /3 = 0.068 ± 0.014. 

We separately cross-correlate each of the three compo¬ 
nents of the sample (SDSS-identified point sources, ex¬ 
tended sources, and non-matches) with the full quasar 
sample and measure a clustering amplitude at 1 h~^ 
Mpc at least 5 a greater than zero for all three cases 
(Figure]^. The clustering amplitude for the SDSS point 
sources and the non-matches are 3 and 4 a lower than 
the clustering amplitudes for the SDSS extended sources, 
respectively, indicating that the non-matches and point 
sources contain more contaminating stars and low red- 
shift galaxies than the extended sources. Nevertheless, 
all three components contain 0.65 < z < 0.9 galaxies, 
and we believe it is prudent to include all three compo¬ 
nents in our sample to eliminate the possibility of seeing- 
dependent variations in galaxy density. 

4. DEPENDENCE OF CLUSTERING AMPLITUDE ON 
QUASAR PROPERTIES 

We measure wlrp) for seven groups spanning 1.3 
decades in luminosity and eight groups spanning 1.3 
decades in virial mass. The seven groups in luminos¬ 
ity consist of three DR7 groups and four DRIO groups, 
while the eight groups in virial mass consist of four DR7 
groups and four DRIO groups (see Figurej^for luminosity 
and virial black hole mass distributions). 

We use linear least squares to find the clustering am¬ 
plitude at 1 h~^ Mpc (/?) and its standard deviation for 
each group, using all five radial bins in each fit. Since we 
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Figure 5. Top left: Angular overdensity w{rp) and two-parameter curvefit w = ^ for the cross-correlation of 4.2 million 

WISE galaxies about 7,049 quasars from DR7 and DRIO. Errorbars are 1 a, and rp is measured in comoving Mpc. Top 
right- Same data as the top left panel, but with the y-axis replaced by rpw(rp) to show that we measure a significant overdensity 
of galaxies even at separations > 1 h~^ Mpc. The curvefit is shown in red. The lower points are from the cross-correlation of 
the WISE galaxies with 37,402 2.2 < « < 3.5 quasars, measured in the same angular bins as the z ~ 0.8 sample. Bottom left-. 

contour for 5 and /3 with 68%, 90%, and 99% confidence intervals. Bottom right: w{rp) divided by the projected correlation 
function of the linear-regime matter field, wp^matter(rp), yielding a quantity proportional to the quasar-galaxy bias &qg- 


do not expect the power-law slope to vary with luminos¬ 
ity or black hole mass, we fix the slope at <5 = 0.84 in each 
fit. If we allow the slope to vary, our results for /3 differ 
by less than 0.3 cr for all luminosity and black hole mass 
subsamples. Table gives the luminosity, virial mass, 
clustering amplitude and for each of these samples, 
Figure l^plots luminosity against clustering amplitude 
for the seven groups in luminosity and virial mass against 
clustering amplitude for the eight groups in virial mass. 

We do not observe any luminosity dependence in the 
quasar-galaxy clustering amplitude. We fit a power law 
between luminosity L and clustering amplitude ff: 

( 12 ) 

The factor in the denominator, 10^^'^^, is the geometric 
mean of the luminosities for each of the seven bins. We 
findp = dlog/?/dlogL = -0.01i0.06, a = 0.105i0.006 
and = 4.75 with 4 degrees of freedom. 

We also find no relationship between quasar-galaxy 
clustering amplitude and virial black hole mass (Fig¬ 


ure Q . We fit a power law between virial mass M and 
clustering amplitude j3: 

( M Y 

= ( 13 ) 

As before, the factor in the denominator is the geometric 
mean of the virial masses of the eight bins. Minimizing 
yields p = 0.02±0.06, a = 0.106±0.006 and = 3.88 
with 5 degrees of freedom. 

Since the continuum luminosity is used to estimate 
virial black hole mass (Equation!^, i-band luminosity 
and virial black hole mass are corrSated (correlation co¬ 
efficient p — 0.40), and the constraints on luminosity de¬ 
pendent clustering and black hole mass dependent clus¬ 
tering are not independent. In order to provide an inde¬ 
pendent constraint on clustering, we also measure the de¬ 
pendence of clustering on Mgll FWHM. If the luminosity 
(and therefore the broad-line region radius) is held con¬ 
stant, all of the variation in line width can be attributed 
to variation in black hole mass. We find that dividing the 
DR7 quasars into four equally sized groups based on Mgll 
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Figure 7. Luminosity dependence {top left), virial black hole mass dependence {top right), and Mgll FWHM dependence 
{bottom) of 13, the clustering amplitude at 1 h~^ Mpc. Luminosity dependence and virial mass dependence are presented for 
both DR7 (black circles) and DRIO (blue squares) together, and FWHM dependence Is presented separately for DR7 {bottom 
left) and DRIO {bottom right). jS is determined from a linear least squares power-law ht to ui(rp) measured in the five radial 
bins of Figurej^with <5 fixed at 0.84. Error bars are 1 a. The best-fit power law is shown in red. The red dashed lines are power 
laws with the same amplitude as the best-fit power law, but slopes ±lcr from the best fits. 



Table 1 


Quasar Sample 

Galaxy Sample 

l^QSO 

logL/L© 

logMeu/M© 

/3±la 


DR7 and DRIO, 0.65 < z < 0.9 

All 

7049 

11.28 

8.49 

0.105 

± 

0.006 

10.23 

DR7 and DRIO, 0.65 < z < 0.9 

Point sources 

7049 

11.28 

8.49 

0.068 

± 

0.014 

2.45 

DR7 and DRIO, 0.65 < z < 0.9 

Extendeci sources 

7049 

11.28 

8.49 

0.133 

± 

0.009 

4.72 

DR7 and DRIO, 0.65 < z < 0.9 

Non-matches 

7049 

11.28 

8.49 

0.098 

± 

0.008 

8.63 

DRIO, 0.65 < z < 0.9, logL < 10.81 logL© 

All 

1052 

10.68 

8.10 

0.119 

± 

0.016 

6.11 

DRIO, 0.65 < z < 0.9, 10.81 logL© < logL < 10.98 logL© 

All 

1051 

10.90 

8.29 

0.097 

± 

0.016 

14.76 

DRIO, 0.65 < z < 0.9, 10.98 logL© < logL < 11.15 logL© 

All 

1052 

11.06 

8.41 

0.100 

± 

0.013 

1.33 

DRIO, 0.65 < z < 0.9, 11.15 logL© < logL < 11.43 logL© 

All 

865 

11.26 

8.54 

0.090 

± 

0.017 

4.49 

DR7, 0.65 < z < 0.9, logL < 11.66 logL© 

All 

1422 

11.55 

8.64 

0.116 

± 

0.012 

1.16 

DR7, 0.65 < z < 0.9, 11.66 logL© < logL < 11.83 logL© 

All 

711 

11.74 

8.71 

0.119 

± 

0.016 

5.90 

DR7, 0.65 < z < 0.9, logL > 11.83 logL© 

All 

710 

12.02 

8.92 

0.084 

± 

0.020 

5.72 

DRIO, 0.65 < z < 0.9, logMen < 8.04 logM© 

All 

1047 

10.88 

7.79 

0.114 

± 

0.016 

3.92 

DRIO, 0.65 < z < 0.9, 8.04 logM© < logALeu < 8.35 log M© 

All 

1047 

10.96 

8.20 

0.088 

± 

0.014 

3.87 

DRIO, 0.65 < z < 0.9, 8.35 logM^ < logMnH < 8.64 log Mf;i 

All 

1047 

11.01 

8.49 

0.106 

± 

0.015 

2.21 

DRIO, 0.65 < z < 0.9, logMen > 8.64logM© 

All 

1048 

11.11 

8.90 

0.100 

± 

0.015 

2.00 

DR7, 0.65 < z < 0.9, logMeu < 8.51 logM© 

All 

703 

11.62 

8.25 

0.106 

± 

0.019 

4.05 

DR7, 0.65 < z < 0.9, 8.51 logM© < logMau < 8.76logM© 

All 

703 

11.68 

8.63 

0.126 

± 

0.018 

2.12 

DR7, 0.65 < z < 0.9, 8.76 logM© < logMeu < 8.99logM© 

All 

703 

11.71 

8.86 

0.098 

± 

0.018 

1.42 

DR7, 0.65 < z < 0.9, logMeu > 8.99logM© 

All 

704 

11.80 

9.20 

0.123 

± 

0.020 

8.74 


Note. — Clustering amplitude (/3), standard deviation, sample size, and mean luminosity and virial black hole mass for the quasar 
subsamples. The mean luminosity and virial mass are weighted using the DRIO redshift distribution (z = 0.77) divided by the redshift 
distribution of each quasar subsample. The clustering amplitudes are derived from the measurements of w{rp) using two-parameter 
minimization for the fit to the combined data (first row) and linear least squares assuming <5 = 0.84 for all other subsamples. 
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FWHM produces four groups with very similar luminos¬ 
ity distributions; similarly, dividing the DRIO quasars 
into four equally sized Mgll FWHM groups produces 
four very similar luminosity distributions. Since the DR7 
and DRIO quasars have very different luminosity distri¬ 
butions, we compute the Mgll FWHM dependence of 
clustering separately for DR7 and DRIO. We also re¬ 
move quasars with Mgll rest-frame equivalent width < 
10 A from our sample (6 in DR7 and 13 in DRIO) be¬ 


cause 


Shen et al. 


is low, the hWH) 


find that if the equivalent width 
ot low S/N quasars is biased low (see 


their Figure 7). For DR7, we find p = 0.24 ± 0.22, and 
for DRIO, we find p = 0.39 ± 0.19 for a power-law fit 
between clustering amplitude and Mgll FWHM. While 
these results hint at a FWHM-dependent trend, more 
data is needed to determine whether clustering is better 
correlated with FWHM than with virial mass. 

Simulations indicate that the luminosity and black- 
hole mass dependence of clustering is st ronger for the 


one-halo term than for the two-halo term (Kauffmann & 
Haehnelt 2002 Thacker et al.|[2009 |. Since we interpret 


the sharp increase in bias at Vp = 0.3 h~^ Mpc (Fig¬ 
ure]^ as indicative of one-halo clustering in this bin, we 
measure the luminosity and black hole mass dependence 
of clustering using this bin only to isolate the luminos¬ 
ity and black hole mass dependence of one-halo cluster¬ 
ing. We do not detect luminosity or virial mass depen¬ 
dence of clustering in the smallest bin: for the luminos¬ 
ity dependence we find p = 0.03 ± 0.11, a = 0.33 ± 0.03 
and = 7.16 with 4 degrees of freedom, whereas for 
the virial mass dependence we find p = 0.14 ± 0.11, 
a = 0.32 ± 0.03 and = 8.54 with 5 degrees of free¬ 
dom. 

To verify the robustness of our results, we repeat the 
measurements of luminosity and virial mass dependent 
clustering using different samples of quasars and galax¬ 
ies. We restrict the galaxy sample first by excluding 
SDSS-identified point sources and then by also excluding 
non-matches. When we exclude point sources, we find 
identical slopes in the /3-luminosity and /3-virial mass re¬ 
lationships. When we exclude both point sources and 
non-matches, we find that the slopes of the /3-luminosity 
and /3-virial mass relations are lower by about 2 a. This 
suggests that our results are not affected by systematic 
errors arising from the sample composition. 

We also study the black hole mass dependence of clus¬ 
tering using virial masses from H/3 instead Mgll. We use 
a restricted redshift range, 0.65 < z < 0.85, to ensure 
that the entire H/3 profile is measured by the DR7 spec¬ 
trograph. We obtain a slope of p = 0.07 ± 0.07 for the 
H/3 masses compared to p = 0.00 ± 0.06 for the Mgll 
masses in this range. Thus, we detect no dependence of 
clustering amplitude on either H/3 or Mgll based virial 
masses. _ 

As discussed in Section [T^ the virial black hole masses 
possess sign ificant uncertainties of up to 0.5 decades 
(Shen|2013|). As a result, the mean true mass of each bin 
is less extreme than the mean virial mass, suppressing 
the virial mass dependence of quasar-galaxy clustering. 
We estimate this suppression by simulating our data us¬ 
ing an assumed power-law relationship between /3 and 
true mass, scattering true mass to virial mass, binning 
the simulated data by virial mass, and fitting a power 


law to the binned data. We then find the ratio between 
the assumed power law slope and the fitted slope. To 
find the distribution of true black hole masses, we ap¬ 
proximate our observed virial mass distribution as log¬ 
normal with scatter CTobs (0.38 decades for DR7 and 0.43 
decades for DRIO) and assume that the distribution of 
virial masses at fixed true mass is lognormal with scat¬ 
ter CTvir; then the true mass distribution is lognormal 
with scatter dtrue = \/o'obs ~ '^vir ■ Previous measure¬ 
ments of the black hole mass function, using virial masses 
but explicitly modelling both the luminosity-dependent 
bias and the Malmquist bias, found atr up = 0-26 and 
gyir = 0.26 for DR7 quasars at z = 0.8 (Shen & Kelly 


20121 . 

We compute the ratio between the assumed power-law 
slope and the measured power-law slope for different val¬ 
ues of CTvir, finding that the ratio reaches 2 at cTyir = 0.33 
(cTtrue = 0.19 for DR7, 0.28 for DRIO) and asymptotes 
as CTvir approaches (Jobs ■ We also note that if luminosity- 
dependent bias is present, atrue > \/^obs ~ ^vir> since 
the luminosity-dependent bias wil l decrease th e obs erved 
scatter in the v irial masses. Since Shen (20131 and Shenf 
& Kelly (2012) find evidence for luminosity-dependent 
bias tor Algll-based masses at z = 0.8, CTyir = 0.3 — 0.35 is 
most consistent with previous measurements of the black 
hole mass function and estimates of Uvir. This value of 
CTvir corresponds to a suppression of the power-law slope 
by a factor of 1.5-2. 

5. DISCUSSION 

5.1. Slope of the eorrelation function 

Comparing the clustering amplitude /3 to previous 
measurements of the correlation length rg is beyond the 
scope of this paper, since it would require converting 
our measurements of the angular overdensity to mea¬ 
surements of the projected cross-correlation function us¬ 
ing the three-dimensional density of WISE galaxies po 
(Equation [^. However, because of the scale indepen¬ 
dence of the power law fit, we can compare our power-law 
slope, 5 = 0.84 ± 0.05, to other results, using Equation |4 
to rela te 5 to other measurements of 7 = d-f 1. [Coil et ah 
(2007) finds 7 between 1.6 a nd 2 depending on the sam¬ 
ple, with 1 CT errorbars ~ 0.2. Porciani & Norberg (2006 1 


find 7 = 1.7^[( p 7 for quasars and galaxies with redshift 
0.8 < z < 1.06’ IShen et al.| (| M1 finds q = 1.69 ± 0.07 
for redshifts 0.3 < z < 0.9, and" Zhang et al. (2013) finds 
7 = 2 . 1 ± ^ 0.2 for redshifts 0.6 < z < 1 . 2 . 

Our measurement of <5 is similar to most of these re¬ 
sults, although very slightly higher. A slightly steeper 
correlation function is to be expected for our sample: our 
color cut biases the galaxy sample in favor of red galaxies, 
which cluster more strongly than blue galaxies. Accord- 
in g to the galaxy auto correlation function measurements 
of Zehavi et al. (2011), red galaxies have 7 = 1.94±0.03, 
whereas blue galaxies have 7 = 1.66±0.03. Furthermore, 
we may measure a steeper slope because we measure 
quasar-galaxy clustering at smaller scales than previous 
studi es. While we mea sure clustering from 0.2-6.4 h~^ 
Mpc, Coil et al. (2007 | measure clu stering from 0.1-10 
h~^ Mpc, Porciani fc Norberg ( 2006|) mea sure clustering 
from 3-20 Mpc, and Shen et al. ^2013 ) measure clus¬ 
tering from 2-25 h~^ Mpc. Therefore, our measurements 
are more sensitive to the one-halo term than previous 
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measurements, leading to a steeper power law slope at 


scales < 1 h ^ Mpc (e.g. Siren et al.||2013[ ). Indeed, our 
results suggest a break in the power law at 1 Mpc, 
with a steeper slope at smaller scales and a shallower 
slope at larger scales (see Figure [5 1. However, we do 
not have enough bins to fit either aiull halo occupation 
distribution (HOD) function or a broken power law. 

5.2. Luminosity dependent clustering: comparison to 
previous results 

We find no dependence of clustering amplitude upon 
luminosity (Figure |7]). We measure the luminosity de¬ 
pendence of clustering across a larger luminosity range 
and with more bins in luminosity than previous studies 
at similar redshifts. As a result, our power-law fit of 
clustering amplitude versus luminosity provides a simi¬ 
lar or tighter constraint on the luminosity dependence of 
quasar clustering than previous results. 

Most previous studies of luminosity dependent cluster¬ 
ing measure the correlation length rg as a function of 
luminosity. Since rg oc (Equation]^, we 

fit a power law between luminosity and to compare 
our results to previous results, using errors computed 
from Taylor series error propagation. We find a slope of 
0.002 ±0.027. A power law is advantageous because it is 
scale independent, allowing us to compare our slope to 
the slopes obtained from power-law fits between rg and 
luminosity. When fitting power laws to other results, 
we assume uncorrelated errors and use the same non¬ 
linear minimization used to determine the amplitude 
and slope of the power law fit to our data. We compare 
to measurements of rg rather than to measurements of 
the quasar bias bg because measurements of the quasar 
bias incorporate errors arising from measurement of the 
g alaxy bias. 


Zhang et al. (2013), measuring the cross-correlation of 
SDSS Stripe 82 galaxies and SDSS quasars at z ~ 0.8, 
found no significant differences in clustering between two 
bins in lumino sity with 10-25% e rror on the clustering 
measurements. Zhang et al. (2013) compared the cluster¬ 
ing of faint and bright quasars tor three different redshift 
bins, yielding six bins in luminosity but with differing 
redshift distributions. In order to fi t a power law to th e 
correlation length measurements of Zhang et al. (2013), 
we normalize their correlation length measurenients in 
the 0.8 < z < 1.0 and 1.0 < z < 1.2 bins so that the 
average correlation length in each redshift bin is identi- 
cal. We obtain a p ower law slope of —0.01 ± 0.10 for the 
Zhang et al. (2013) measurements, for quasars at z = 0.9 
a nd 10.98 < iog±(g < 11.57. 


Shanks et al. (2011) compared correlation lengths mea- 
sured by SDSS, 2tjZ, and 2SLAQ at z = 1.4 across 
slightly more than a decade in luminosity (11.44 < 
logL© < 12.48). A power law fit to their results yields a 
slope of 0.04 ±0.03, providing a similar constraint on the 
magnitude of luminosity-dependent clustering as our re- 
sult, although the measurements of [Shanks et al. (2011) 


a re at higher r edshif t than our resurtsT 
jShen et al. (2013), measuring the cross-correlation of 
DRIO CiVlASS galaxies and DR7 quasars at 0.3 < z < 0.9 
(z = 0.53) with 11.19 < log^Q < 11.88, found no lumi¬ 
nosity dependenc e over four bins in luminosity. We fit a 
power law to the Shen et al. (2013) measurements of lu¬ 


minosity and correlation length, obtaining an exponent 
p = —0.01 ± 0.05. While the cl ustering measurements 
in each bin of Shen et al. (2013) have « 10% errors in 
^(r), compared to our errors of « 15% in w(rp), we pro¬ 
vide a tighter constraint on the luminosity dependence of 
clustering because we measure clustering across a larger 
luminosity range and with more bins in luminosity. 

Comparison of our work to previous work shows that 
measuring clustering across a wide range of luminosity 
is important in obtaining a more precise measurement of 
(or constraint on) luminosity-dependent clustering. The 
restricted range in optical luminosities is often cited as 
a justification for the lack of mea sured luminosity de¬ 
pendence (e.g. Bonoli et al. 2009) and the wider lumi¬ 
nosity range probed by A-ray selected AGN may explain 
why X-ray luminosity dependent clustering has been de- 
tect ed while optical lumi nosity dependent clustering has 
not ( Krumpe et al.|[2012 ). 


5.3. Luminosity dependent clustering: comparison to 

models 

In Figure [^ we compare our results to theoretical pre¬ 
dictions of luminosity-dependent cl ustering from four dif- 

(|2009|), referred 


ferent mo dels of quasar evolutio n: Shen 


to as S 09, Hopkins et al. (2014| (H14), Conroy & White 
(|2013|) (CW13) , and a simple “lightbulb” model used by 
Hopkins et al. ([2007 ) for comparison to more sophisti- 
T 


cated models, 'i'he S09 and H14 models are physically 
motivated, while the CW13 and light bulb models merely 
provide a relationship between halo mass and quasar lu¬ 
minosity. 

The “lightbulb” model is the simplest of the four con- 
side red. We use the relat ionship given at z = 1 in Figure 
3 in Hopkins et al. (2007) (the line labelled “All QSO at L 
= LEdd”)- In this model, all quasars radiate at the same 
Eddington ratio and black hole mass is related to host 
halo mass without scatter. Thus, the luminosity-bias re¬ 
lationship is entirely determined by the halo mass-bias 
relationship. However, observation of scatter in both the 


et al. 2006|) and 

in the halo mass-black hole mass relation- 

ship (I'errarese 

2002 

) indicate that the lightbulb model 

is an oversimp 

litication. 'i'he CW13 model (using the 


z = 1 relationship from Figure 8) is a “scattered light- 
bulb” model with scatter in both the Eddington ratio 
distribution and in the black hole mass-host halo mass re¬ 
lation. Similar to the lightbulb model, the CW13 model 
uses a step function for the quasar light curve and uses 
the quasar lifetime to set the quasar bias at a given red¬ 
shift . 

In contrast, both the S09 and H14 models are based on 
physically motivated light curves. The H14 model (using 
the z = 1.5 bias-luminosity relationship in their Figure 
5) considers two modes of quasar fueling: major merger 
induced fueling, which is responsible for most high lumi¬ 
nosity quasars, and stochastic fueling mechanisms snch 
as minor mergers, which are responsible for lower lu¬ 
minosity quasa rs. The merger-induc ed fueling model is 


quite similar to Hopkins et al. (2007), and predicts from 
merger simulations that host halo mass and p eak quasar 


luminosity are tightly linked ( Lidz et al.|2006 ). However, 


the relationship between peak luminosity and instanta¬ 
neous luminosity has considerable scatter, arising from 


the merger-based quasar light curve of Hopkins et al. 
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(20051. H14 also considers stochastic accretion from low- 
luminosity quasars and Seyfert galaxies, which reside in 
disk hosts rather than bulges, although the fraction of 
disk-hosted quasars at the luminosities of our sample is 
estimated to be quite small. In the H14 model, quasars 
spend most of their lives accreting at low rates. As a 
result, many faint quasars have high peak luminosities 
and thus high host halo masses. Thus, a large range in 
instantaneous luminosity maps to a small range in peak 
luminosity (and thus black hole mass, halo mass, and 
clustering strength). 

The S09 model (using the relationship at z = 1 from 
their Figure 5) is similar to H14, differing in only three 
respects. First, S09 estimates the major merger rate us¬ 
ing simulations while H14 estimates the merger rate from 
empirical halo distribution models. Second, H14 derive 
their quasar light curve from simulations of galaxy merg¬ 
ers, while S09 fits the light curve empirically from obser¬ 
vations. Last, S09 does not consider alternative fueling 
mechanisms for low-luminosity quasars, but, as discussed 
in the paragraph above, the practical significance of this 
choice is small. 

Each model specifies the relationship between the lin¬ 
ear quasar bias bq and the quasar luminosity. In analogy 
with Equation the linear bias bq is defined using the 
quasar-galaxy cross-correlation function: 



Figure 8. Measured clustering amplitude as a function of 
luminosity, compared to theoretical predictions. All four the¬ 
oretical models predicted quasar bias as a function of luminos¬ 
ity. We assumed that quasar bias is proportional to clustering 
amplitude and normalized the relationship to the correlation 
amplitude of the median luminosity bin. Luminosities are 
measured using Mi{z — 2), measured in solar luminosities, 
corresponding to rest-frame 2500 A luminosity. 


bq = 


^QG 


bc^n 


(14) 


If a linear quasar bias provides an appropriate fit to our 
data, then Equations m and imply that fi o^bq. To 
compare our results to models predicting the luminosity 
dependence of quasar bias, we convert bq to wirp) by 
forcing each model to match the median /3 at the median 
luminosity of our data. 

The linear quasar bias is only applicable on larger 
scales. At scales < 1 h~^ Mpc, the one-halo term leads 
to nonlinear clusteriM and deviations from Equation |14[ 
as shown in Figure However, Figure implies that 
a linear bias is appropriate in all bins except the inner¬ 
most bin, implying that P oc bq is not an unreasonable 
assumption. Moreover, we expect that one-halo cluster¬ 
ing will also depend on lu minosity, perhaps more strongly 


than two-halo clust ering (Kauffmann & Haehnelt 2002 
Thacker et al.|[2009 l. 


i'he models present luminosity in different units. The 
S09 model gives the relationship between bias and bolo- 
metric luminosity in erg s~^. We conv ert to Mi(z = 2) 


using Equation 1 from Shen et al. (20091 and then to solar 
luminosity using the same methods as in this paper. The 
H14 and lightbulb models express luminosity in terms of 
solar bolometric luminosity; we convert to erg s“^ using 
the bolometric luminosity of the sun, Lq = 3.827 X 10^' 


c-l 


et al.| (20091 to convert to Mi(z = 2) and the same rela- 
in: 


As before, we then use Equation 1 from Shen 


tionship used in this paper to convert to solar luminos¬ 
ity. The CW13 model gives luminosity using Mj(z = 2), 
which we again convert to solar luminosity following the 
methods used above. 

The H14 model provides the best fit to our data (x^ = 
5.96 for 6 degrees of freedom), followed by the S09 model 
(x^ = 6.43), the CW13 model (y^ = 10.38), and the 
lightbulb model (x^ = 11.22) (Figure^. Our data prefer 


the H14 model to the lightbulb model at 2 cr. However, 
at present none of these models are disfavored by our 
data at more than 2 a. The divergence between models 
is greatest at high luminosities, suggesting that improved 
constraints could be obtained by measuring clustering at 
higher redshifts, or by using a deeper sample of tracer 
objects to obtain more precise clustering measurements. 

5.4. Black hole mass dependent clustering 

We find no significant relationship between virial black 
hole mass and clustering strength, in contrast to previous 
measurements of blac k hole mass dependen t clustering by 
Fine et al. (2006) and Komiya et al. (20 13|). Our method s 
are quite difterent from the methods mFme et al. (20061; 
they measured the relationship between black hole mass 
and host halo mass for bins at different redshift, allow¬ 
ing them to pro be a larger dynam ic range in black hole 
mass. Moreover, Fine et al.|(2006 1 measured only a weak 
relationship between black hole mass and host halo mass. 
Our results are in mo re significant tension with those of 
Komiya et al. (20131: a power-law fit to the black hole 
mass d ependence of the co rrelation length ro from the 
data of Komiya et al. (20131 yields a slope of 0.14 ± 0.05, 
compared to a slope of 0.01 ± 0.03 for our data (fitting 
the black hole mass dependence of cx tq). There are 
significant m ethodological differen c es between our results 
and those of Komiya et al. 1 (|2013||: |Komiya et al.| (|2013|) 
measured black hole mass over a signihcantly broader 
redshift range (0.1 < z < 1.0), used masses from three 
virial mass indicators (Mgll, H,5, and Ha), and did not 
correct for differing redshift distributions between differ¬ 
ent bins in black hole mass, although they found only 
minor differences using reselected bins containing identi¬ 
cal redshift distributions. We also note that other studies 
failed to find a significant relati onship between bla ck hole 
mass and clustering strength: Shen et al. (|2009 1 found 
no dependence of quasar clustering on virial mass using 
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two bi ns in viral mass at 0.4 < z < 2.5, and Zhang et al.| 
(20131 found a 1-2 cr difference in clustering strength 
between two bins in virial mass at 0.6 < z < 1.2. Fur¬ 


thermore, Chatterjee et al. (20131 found no significant 
relationship between host cluster mass and virial mass 
at 0.1 < z < 0.3; rather than using clustering to mea¬ 
sure the linear bias, they directly matched quasars with 
galaxy clusters and estimated the cluster mass using the 
optical richness. 

Our results also appear to conflict with the predic¬ 
tions of many different models, which predict a rela- 


(ILidz et al. 2006 

Kauffmann & Haehnelt 

20021 

& White 


h’errarese|2002|. These moc 

els all 


tiic i cLii^c o»i iiictooco oo»iioioici coi iicic, 

with 1.2 < a < 1.8. We compare our data to these pre¬ 
dictions by estimating the ,d-black hole mass relationship 
that would result from these Mbh~ A/dmh relationships. 

We begin by relating Mdmh to the linear quasar b ias 
bq. We use the fitting formula of Sheth et al. (20011 to 
find 6 (M): - - 


6=1 


y/a5sc[z) 


\/a{ai'^) + \/ab{ai'^Y 

_ {av^Y _ 

{av'^Y + ^(1 ~ c)(l — c/ 2 ) 


(15) 


where a = 0.707, 6 = 0.5, c = 0. 6, v = 6sc{z)/(t(,M, z ) 
and Ssciz) = 0 . 15 ( 127 r) 2 / 30 )/o °55 (|Navarro et al.||l997|. 
a{M,z) = a{M,z = 0)D{z), where L>[z) is the linear 
growth factor. We approxi mate D{z) using the fitting 
form of Carroll et al. ( i992[ ): 


D(z) = 
Di{z) = 


Di{z) 


1 + 


Di(0) 

1 

iTz 




(16) 


where flmz and Daz give the evolution of the cosmolog¬ 
ical parameters and IIa. ct(M, z = 0 ) is the fluctua¬ 
tion in the density field for a halo of mass M at redshift 
0 , given by 

1 

/ k'^p{k)W'^{kr)dk (17) 
ZTT Jo 

where P{k) is the power spectrum of the linear-regime 
matter field, computed from the transfer function of 
Eisenstein & Hu ( |1998[ ), 


W(kr) = 


3 [sin kr — kr cos kr] 

{krY 


(18) 


and r{M) = (3M/(47rpo))^'^^, where po is the mean den¬ 
sity of the Universe at z = 0, 2.78 x 
Mpc“3. We normalize the linear power spectrum such 
that cr(r = 8h~^Mpc) = 0.84. 

Since our angular clustering measurement does not al¬ 
low for accurate estimates of bq, we relate bq t o P us- 
ing previous measurements of bq at z = 0.8. Groom 


et al. (2005) measured bq{z = 0.8) = 1.49 ± 0.21; they 
also provided an empirical quadratic fit to their mea¬ 
surements of bq from z = 0.53 to z = 2.48, yielding 
bq{z = 0.8) = 1.47 ± 0.2. We fi t a quadrati c function to 


the measurements of Shen et al. |2009 201^ White et al. 


(|2012|,[Padmanabhan et al .|(|2009t 
Forciani lNorbergl2006j, using 


Groom et al.| ( |200^ 


9 different nieasure- 
ments with 0.3 < z < 3.8 and obtain bq{z = 0.8) = 1.43. 
Following these estimates, we use bq = 1.5 as the linear 
bias of our sample. Since bq oc /3, we multiply b{M) by 
the ratio between /3 = 0.105 and b(z = 0.8) = 1.5. Since 
b{M) is approximately linear over 1 decade in host halo 
mass, /3(Mbh) is well-approximated by a power law. 

The power-law slopes 1.2 < d log Mbh/cI log Mdmh < 
1.8 correspond to 0.1 < d log/3/d log Mbh < 0.15. We 
only find « Icr difference between these slopes and our 
measurement of the Mgll virial mass dependence of clus¬ 
tering, after adjusting for the factor of 1.5-2 suppression 
arising from the large scatter in the virial masses. Our 
results are also consistent with no relation ship between 
Mbh a nd Mdmh, a possibility suggested by [Shanks et al.| 
(2011 1 , who argue that the lack of luminosity-dependent 
clustering and the corresponding long quasar lifetime 
supports a pure luminosity evolution model, in which 
bright quasars at high redshift become faint Seyfert 
galaxies at z = 0. In the pure luminosity evolution 
model, host halo mass is not tightly rel ated to black hole 
mass, as s upported by the results of [Kormendy et al(1 
(2011) and Kormendy & Bender] (2011), who found no re- 
lationship between black hole mass and stellar disk mass 
in disk galaxies at z = 0. Our results suggest that dis¬ 
criminating between models with no link between black 
hole mass and host halo mass and models with a tight 
relationship between black hole mass and host halo mass 
will require clustering measurements with larger samples 
and larger dynamic range in black hole mass. 


6. CONCLUSIONS 

Measuring the luminosity dependence of the quasar 
clustering amplitude allows us to test different rela¬ 
tionships between host halo mass and quasar luminos¬ 
ity. Previous measurements of the quasar autocorrela¬ 
tion function and the three-dimensional quasar-galaxy 
cross-correlation function have suffered from small sam¬ 
ple sizes: the spatial density of quasars is low, and it is 
difficult to obtain a large number of spectroscopic red- 
shifts distributed across the sky. This study alleviates 
these concerns by measuring the angular overdensity be¬ 
tween galaxies and quasars, resulting in a much higher 
density of the tracer population. 

We find no luminosity dependence of the quasar-galaxy 
cross-correlation function, consistent with previous find¬ 
ings. A power law fit of luminosity to clustering am¬ 
plitudes gives a slope of — 0.01 ± 0.06, a much tighter 
constraint than those provided by previous studies at 
similar redshift. We also fail to detect a relationship 
between clustering strength and Mgll-based virial black 
hole mass. However, this result is consistent with the¬ 
oretical predictions tightly linking black hole mass and 
host halo mass. These results indicate that, within the 
ranges of luminosities considered here, the most luminous 
quasars reside in a wide range of dark matter halos. 
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